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Abstract 

We study cosmology in the bigravity formulation of the dRGT model where matter couples to both 
metrics. At linear order in perturbation theory two mass scales emerge: an hard one from the dRGT 
potential, and an environmental dependent one from the coupling of bigravity with matter. At early 
time, the dynamics is dictated by the second mass scale which is of order of the Hubble scale. The set 
of gauge invariant perturbations that couples to matter follow closely the same behaviour as in GR. 
The remaining perturbations show no issue in the scalar sector, while problems arise in the tensor 
and vector sectors. During radiation domination, a tensor mode grows power-like at super-horizon 
scales. More dangerously, the only propagating vector mode features an exponential instability on 
sub-horizon scales. We discuss the consequences of such instabilities and speculate on possible ways 
to deal with them. 


1 Introduction and Summary 

The physical mechanism responsible for the present day acceleration of our universe is unknown. The 
simplest explanation is a positive cosmological constant; however the large amount of tuning required to 
fit the data seems excessive. Alternatives based on modifications of Einstein General Relativity (GR), 
which become observationally relevant at large scales, are being actively explored nowadays [T]. Among 
them. Massive Gravity [2] has received special attention: from an effective held theory perspective, it is 
one of the most natural options to investigate when renouncing to the diffeomorphism invariance of GR. 
In such a scenario, the graviton mass introduces a new energy scale that can be related with the scale of 
dark energy. In its simplest incarnation, to build a massive deformation of GR a reference non-dynamical 
metric is needed. Besides the aether-like nature of the reference metric, an unattractive feature from a 
theoretical perspective, there are various motivations to go beyond massive gravity and enter in the realm 
of bigravity theories [aHElE]. For example, spatially hat homogenous Friedmann-Robertson-Walker 
(FRW) solutions do not exist [3 in Lorentz invariant ghost free massive gravity, and even allowing for 
open FRW solutions [8] strong coupling [9] and ghostlike instabilities [TO] develop. Flat FRW solutions 
exist mm in the case of Lorentz breaking models (nidiKEKiniiiTi. 


In bigravity one can find various branches of regular cosmological solutions describing flat FRW 
cosmologies [lanH]. A branch where the gravity modification is equivalent to an effective cosmological 
constant suffers of strong coupling [20]. A more promising branch is unstable at early time, being 
characterized by an exponential growth of fluctuations Also, singular FRW-type solutions 

exist [T8| which exhibit only mild instabilities (power-law growth of vector and tensor modes [221123j l. 
On the other hand, they correspond to bouncing universes characterized by a naked curvature singularity, 
which makes their physical relevance questionable. Recently, it has been proposed to extend the theory 
of massive (bi) gravity by considering a more general coupling to matter - called doubly matter coupling 
- in which the physical metric coupled to the matter energy momentum tensor is an appropriate linear 
combination of the two metrics |24[l25j (see also |26II27) for different approaches). This scenario, although 
problematic for massive gravity [28] , is potentially interesting in the bigravity setup since a qualitatively 
new branch of FRW cosmological solutions exists [29], hence its cosmological perturbations deserve to 
be investigated. Such a theory is not ghost-free [301 ED, but there exist physically interesting situations 
where the Boulware-Deser (BD) ghost does not represent an immediate phenomenological problem. This 
is the case if its mass is above the cut-off scale Ac of the theory under consideration. In addition, Ac might 
be parametrically larger than the strong coupling scale where the effects of the graviton mass term become 
important and interesting. Moreover, the ghost does not manifest itself at linear order in an expansion 
in fluctuations around particularly symmetric configurations (as for example FRW cosmologies). In this 
work, after studying the two branches of cosmological solutions at the homogeneous level, we focus on 
the dynamics of linearized cosmological fluctuations around the new homogenous backgrounds allowed 
by the doubly matter coupling. No hints of BD ghost mode are found at linear level in fluctuations, 
and the theory propagates the seven degrees of freedom as expected for a healthy bigravity theory. The 
dynamics of scalar fluctuations is healthy, and no instabilities are found in this sector. The dynamics 
of tensor and vector fluctuations is richer, but it shows problematic behaviours. The tensor sector 
exhibits a power-like instability at superhorizon scales during the radiation domination era. We argue 
that such instability is not extremely serious, and can be tamed by an appropriate choice of initial 
conditions, possibly motivated by inflation. Much worse is the behaviour of vector fluctuations. In this 
case, during radiation domination, we find a gradient instability at subhorizon scales, which leads to an 
exponential growth of small scale fluctuations, rapidly driving the theory outside the regime of validity 
of perturbation theory. Hence, this serious instability rules out the cosmological configurations that we 
consider. Nevertheless, we speculate on possible extension of the bigravity theory under consideration, 
that might be able to cure such instability problems. 

2 The theory under consideration 
2.1 Scalar Field 

Before plunging in the study of bigravity, it is interesting to understand in a simplified setting the 
peculiar feature of the non-minimal coupling of matter to gravity proposed by Ref. |24) . We will show 
how the consistency of the effective description of matter as a (perfect) fluid necessarily requires the 
dynamical character of the second metric, selecting bigravity as the only consistent formulation. Take as 
matter a scalar field cj) that couples with gravity not simply by the metric but trough a combination 
of and a non-dynamical flat metric 

ds'f = ffj,udx^dx'^ = —z'‘^ dt^ + dx^dx^Sij = d^^^'diy^^rjab dx^dx^ . ( 1 ) 

We denoted by ' the time derivative with respect to t. In order to restore diffeomorphism (diff) invariance, 
the non-dynamical metric / can be written using four Stuckelberg fields; in the unitary gauge we have 


1 


( 2 ) 


= z{t) and 5*^. The scalar field (f), with a potential F{(j)), couples to gravity according to 

2M2 + . 

Thus, (f> is minimally coupled to G^i, defined by 

+ 2 a /3 , (3) 

where Y'i = {^fX)y and Xy = g^^ ftju and a , (3 are two arbitrary constants. Notice that (j) is not 
minimally coupled to g^y. Setting /3 = 0 and a = 1, we recover the standard minimal coupling. It is 
clear from the ADM canonical analysis that the total Hamiltonian is not anymore linear in the lapse N 
and shift of the dynamical metric g^y, then the number of propagating degrees of freedom (DoF) will 
be more than 3. The seemingly innocent action ([^ actually represents a modification of gravity. 

Let us consider FRW homogeneous cosmological solutions where 

ds^ = dt^ + a^(t) dx^dx^Sij , (4) 



and 


dslg = Gpydx^dx^ = —N^gdt‘^ + a^^dx^dx^5ij ; 
A'eff = a N + 13 z', Oeff = aa + (3 . 


The Energy Momentum Tensor (EMT) for the scalar is diagonal, fluid-like and can be written as 


Tq — Pff, ~ 


T}=p^5), 






2 N 1 , 




a 


XeS CL 


2 

eff 


Na^ 


A '2 


2iVe^ff 


- V 


( 6 ) 


which has a peculiar dependence on N and z'. The expression for the EMT reduces to the standard one 
when /3 —)■ 0, a —)• 1. 

Contrary to the case of a scalar field minimally coupled to gravity, the time-time and the spatial 
components of the Einstein equations and the equation of motion for (p are all independent. Indeed, 
taking the time derivative of the time-time component of the Einstein equations and using the equation 
of motion of (j), one can solve for a"; then inserting this expression in the spatial components of the 
Einstein equations one gets the following constrainlQ 


Pp<t> = 0- (7) 

Thus, unless 13 = 0, Einstein equations require that p^p = 0. The same constraint follows from the 
requirement that the scalar EMT is conserved. Of course such constraint has no counterpart in GR, 
where the EMT for (p is automatically conserved when p satisfies its equation of motion. 

Hence: 

- in a FRW background the dynamics of a scalar field in Q is not equivalent to a perfect fluid; 

- the pressure has to vanish. 

'^We do not consider nnphysical cases where z' and/or aes — 0. 
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Such difficulties were not taken into account in |32j . 

In the following we will show that both issues are absent when the non-dynamical Stuckelberg met¬ 
ric 0 is promoted to a full-fledged dynamical one, see also |28j . The reason for such a behaviour can 
be traced back to the non-dynamical nature of the metric used in the new matter coupling. The 
conservation of the energy momentum tensor 

v(9)r^^ = 0, (8) 

defined as the response of the matter action Smatt to a diffeomorphism variation of the dynamical metric 

<55„,att = - 2 / ’ (9) 

gives such a strong constraint that in general cannot be satisfied unless very special conditions like Q 
are met. This is not very surprising and it is typical of theories with non-dynamical object |33l I34j . 

When instead the metric gets dynamical, inserting its own Ricci scalar in the action, for the 
following FRW parametrization 

ds'j = — dt^ -|- cu^(f) dx^dx^dij , (10) 

condition 0 becomes 

P {z'a'- Nuj') = 0 , (11) 

and a new possibility of realizing 0 opens up. Within this new way, the scalar field dynamics can be 
still captured by the perfect fluid description and no spurious constraint is required. 

Therefore, for the rest of this paper, we will parametrize the matter content of the Universe through 
a perfect fluid and the metric entering in ([^ will become dynamical in the bigravity formulation. 


2.2 Bigravity and Matter Coupling 


Consider the action of massive bigravity with the dRGT potential as interaction between the two dy¬ 
namical metrics and 


S= fd*x{^ [m 2 (7^ - 2m2 U)] + T/kM^ • 


( 12 ) 


In the presence of two metrics, it is not a priori clear to what metric matter couples to. In general, the 
BD ghost revives in the presence of doubly coupled matter [301 ET] . In |21j a new matter coupling was 
proposed where matter is minimally coupled to the effective metric given by Q. Although the BD 
ghost persists even with this special doubly coupled bigravity model, it was shown that the BD ghost 
does not appear in the decoupling limit [241130) . As for the scalar field non-minimally coupled to gravity 


of section 2.1 the theory described by the action (12) propagate more than the seven DoF expected 


in the bigravity formulation of dRGT. The extra scalar mode is most probably a ghost, however, the 
theory is still acceptable if the mass of such a mode is above the ultraviolet cutoff Ac. Clearly this point 
deserves further investigation. 

The matter EMT is defined as the response of the matter action to a variation of g (/) 


5S, 


matt 


= -^j d^x^T^'^dG^, I 


d^x 


fiiy 




fliy 


(13) 
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The modified Einstein equations can be written as 


Ei^ + Q':; = ^Ti:, (14) 

+ = (15) 

where Q (Q) are the effective energy-momentum tensors induced by the interaction term for the two 
metrics. 

Let us introduce the BD ghost free potential ISaiMlEZ! 


4 

^ ^ , 

n=0 


where the Vn are the symmetric polynomials of Y 

Eo = 1, El = ri, E 2 = Ti - r 2 , E 3 = rf - 3 n r 2 + 2 ra , 
E 4 = Ti - 6 rf r 2 + 8 ri Ts + 3 r| - 6 r 4 , 


(16) 


(17) 


with Tn = tr(y"'). We have that 

QO = - (E' E)C] , (18) 

QM = {V Y)>i, (19) 

where (E')0 = dV/dY^ and q = detX = det(/)/det(g'). 


3 Homogeneous cosmological solutions 

3.1 FRW Ansatz and Conservation Laws 

Let us consider homogeneous FRW bi-diagonal metrics with flat spatial slices in conformal time, so the 
form of g and / is as follows: 


ds^ = a^(r) [—dT"^ + dr"^ + 

ds = dr"^ + dr^ + dVl^]^ 

The effective metric gets the following form 

dSgff = —{a a -b l3cio)‘^dT'^ + {aa + /3io)‘^ {dr^ -b dO^) . 
Consistency of the equations of motion requires the following Bianchi-type constraints 

(2go-M-2r,^) =0, 


When T/xu is the EMT of a perfect fluid, i.e. Tfiu = {p + p) U/i + p G^u with Ua up = 
previous equations can be combined to give 


( 20 ) 

( 21 ) 

( 22 ) 
— 1 , the 


3(rc + 1) {a aH + j3ijj T-Lu]) p + {aa + fiu) p' = Q , 

{cH — T-Luj) w afd {aa + ( 3 o;)^ P — 2 m?Mpi (ai -b 4 02 aca -b 6 03 


= 0 . 
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(23) 

(24) 






where T-L = a'/a and = oj'/uj and the equation of state w, p = w p. Notice that condition (23) 


corresponds to the conservation of the matter EMT 'T^'^ with respect to the metric Gn^. The constraint 


(24) can be realised in two inequivalent ways: 


Branch 1 


In this branch eq (24) is realised through the vanishing of the square bracket, i.e. the pressure of the 
fluid is determined by the massive potential. Notice that at early time, when p ;§> consistence 

requires that rc ~ 0, even though the Universe should be radiation dominated at that epoch. In this 


branch we have to deals with the very same issue found in section 2.1 and the description of matter as 
a perfect fluid is not consistent with the one of a scalar field. This branch was studied in [32] , when 
is a flat non-dynamical metric, in presence of a scalar field. Though, contrary to the case of massive 
gravity with minimally coupled matter, flat FRW solutions exist, the non-physical requirement of re = 0 
makes the present branch not very appealing as discussed in Ref. 


Branch 2 


In this case 


c = 


n 


(25) 


Notice that the limit /I —>■ 0 exists and we recover the very same branch of FRW cosmology in bigravity 
with standard matter minimally coupled to the metric Contrary to branch 1, since condition (25) 


is matter independent, the scalar field dynamics is equivalent to a perfect fluid; in this sense matter has 


the standard effective description. R is also interesting to note that (25) is equivalent to requiring that 


the matter’s action and the interaction part of the action are separately diff invariant, namely 

= 0 /is on-shell, 

= 0 / and cj) are on-shell; 


(26) 


and equivalently for the metric /. Clearly, the branch two is the most interesting one and from now on 
we will focus on it. 


3.2 Background solutions for Branch 2 

Introducing the ratio of the two scale factors ^ = uj/a, the tt-component of the modified Einstein 
equations for g reads 




(®o/3 + ai^ + 2a2f + 2a3 f) 

pi 


(27) 


The relation between ^ and a is determined by using (25) and (27) in the tt-component of the modified 


Einstein equations for /; the result is an algebraic equation 

2 m^ [ai + {6 a 2 — K oq) ^ 3 (6 03 — k oi) 6 ( 4 a 4 — k 02) — 6 k 03 

= M~i^p{a + P {ko^ - P) . 


(28) 


The case a —)• 0 and /3 —)• 0 has been already studied in |20| and instabilities were found in the scalar 
sector. 

Throughout this paper we assume that the scale of the graviton mass m is of the order of the present 
Hubble scale, i.e ~ p\ ~ Hq. Such a choice is the most natural one if massive gravity has 
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anything to do with the present acceleration of the Universe. Then, according to the evolution equation 
for the matter energy density (23) that gives 

p = po[ia + /3^)a]-^^^+^'> , (29) 

at early times, provided that ^ <C l/o, we can always consider p <C 1 as a dimensionless expansion 

parameter. Therefore, at early time, solutions of equation (28) can be classified according to the different 
regimes of i.e ^ <C 1, ^ 3> 1, and ^ ~ 1. 

► C < 1 

In this case the solution for the matter energy density is given by 




+ 


2 m? ^ [k ao a /3 — ai (k — 3 /3^) — 6 02 a /?] 


(30) 


a4/32 


+o(a- 


At the leading order we have p ^ p^. This rules out the small ^ regime at early times. 
► ?> 1 

In this case we have that at early times the matter energy density is given by 


p = —l2m^M. 


pi 


a 


12 m? Mpi [«; 02 a /I — 03 (3 K — / 3 ^) — 4 04 a / 3 ] 

K ^ 


(31) 


+ o{^ 


-2\ 


The behaviour is equivalent to a cosmological constant; thus, also the large ^ regime is not suitable 
for early time cosmology. 

► ^ ~ 1 

For this last case we have two solutions for and at the leading order they read 


^={ 


--p+o 
-^ + 0 

K,a 


Pi 


(32) 


When ^ ~ ^ the spatial components of the effective metric (21) are singular and moreover the 

early time cosmology is dominated by a cosmological constant, i.e. oc m?a‘^. When instead 
^ ~ at the leading order we have 


c = 1, 


= 


K a- 


+i}p 








(33) 


and the corrections are of the order ~ 0{m?M‘^Jp). Therefore, up to a renormalisation of the 
Newtonian constanlQ the early time cosmology is very similar to GR at the background level. 

^We should compare the coefficient in front of p in ( |33| with the one appearing in spherically symmetric solutions 
through the Vainshtein mechanism using this new matter coupling. 
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The dynamics of the background is described by the following equations 


din a 
dp 

din a 


(c-l)e, 

-3 (l ++ 


( 34 ) 


where the first equation follows from the definition of ^ and the second equation is the continuity 
equation for matter. The lapse function in the second metric, c, can be found combining the 
Einstein equations and can be expressed in terms of The late time fixed point is given by c = 1 
and p = 0, which corresponds to a de Sitter phase. At early times, ^ ~ ^ and c = 1. Once 
the density becomes lower, the solutions are attracted towards the de Sitter fixed point. This 
transition happens when rin?Mpi/p ~ 1. At the de Sitter fixed point we have ~ this is why, 
in order to explain the late time acceleration of the Universe, we need to assume m Hq. Fig. 1 
shows an example of the two dimensional phase plain spanned by p and ^ and the trajectory of 
the background solution that connects the early time cosmology to the de Sitter fixed point. 



Figure 1: The phase plane described by eq (34). We take oq = ai = 1, 02 = 03 = 04 = 0, a = /3 = 
1 , K = l,r(; = l/3asan illustration. 


Notice that the ^ ~ 1 regime is incompatible with the /3 —)• 0 limitIndeed, when (3 is very 
small, ^ ~ O (jn^Mpi/p'^ and we turn back to the small ^ regime studied in [20]. Thus the early 

time cosmology for ^ ~ ^ exists only when /3 ;§> m?‘M‘^Jp. In order to have this early time 
cosmology until the solution reaches the de Sitter point, (3 needs to be 0(1), assuming that all 
other parameters are also 0 ( 1 ). 

In the next section we will analyse cosmological perturbations around this background. 


4 Cosmological Perturbations 

Perturbations around a FRW background can be classified according to representations of the 50(3) 
group, namely scalar, vector and tensor modes. It is convenient to use the gauge invariant formulation 
following |2ilEI| and for the benefit of the reader the relevant expression are collected in Appendix 

^We stress that in the discussion of the various regimes we have considered /3 7 ^ 0 . 
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In order to avoid cluttering in the main text, the general expressions of the perturbed Einstein 
equations for the scalar and vector sector are given in Appendices B.l and B.2, see also [39] for a 
different derivation. 

The general feature that emerges is that comparing with the ghost free bigravity theory with mini¬ 
mally coupled matter [20|, the new matter coupling gives rise to an additional environmental dependent 
mass', namely, the new equations can be obtained from the old ones by the replacement 


m 


J(a , c, 




(35) 


where J is a function of time that depends on the field under consideration. The bottom line is that two 
mass scales are present; a hard mass coming from the interaction potential for the metrics, see (12); 


and, as a consequence of the non-minimal coupling, a second “soft” running mass scale. In the physical 
branch, where ^ ~ 1, at early times the soft mass always dominates on the hard one, since p/M^^ 3> m?. 
Thus at early time the potential V plays no role and can be simply neglected at the leading order in 
~ O ^m^Mpi/p^. Of course this is the case only when f3 ^ 0, otherwise the regime ^ ~ 1 simply does 

not exist and the corrections become the leading part. The perturbed equations of motion can be 
fully solved in this approximation. 

Notice that the new matter coupling is driven by the matter EMT , so the contribution to 
the ij component of the perturbed Einstein equations is always proportional to w, as a result such a 
contribution vanishes when the Universe is dominated by non-relativistic matter where tc = 0. 

In the following we will often use the Fourier transform with respect to x* of the various perturbations, 
the corresponding 3-momentum will be denoted /c* and k'^ = k'^k^dij. To keep notations as simple as 
possible we will use the same name for the field and its Fourier transform. 


4.1 Scalar Sector 


In the scalar sector, see Appendix |B.1[ the fields S, Bi and 'I'i /2 are non dynamical and can be expressed 
in terms of ^ij 2 ^ basically one of the Bardeen potentials for g and /, that satisfy two second order 
equations; thus 2 scalar DoF propagate. As already pointed out in section |2T] by using canonical analysis 
an additional scalar is expected, however as matter of fact such a mode does not propagate around a 
homogeneous background. Nevertheless it is expected to appear in a less symmetric background and/or 
at higher order in perturbation theory m- The consequences of the missing BD mode, at leading order 
in cosmological perturbations around FRW cosmologies, deserve further study. 

In order to solve for the two propagating fields, it is convenient to use the following combination of 
themselves: 

/32 ^ ^ f 3 ‘^ 


^+ = ^1 + 


K 






K 


r^2 • 


(36) 


It is <I>_|_ that couples to matter thus <I>+ is the relevant metric perturbation for observations. turns 
out to decouple completely (at the leading order) and satisfies the same equation as in GR. Instead the 
equation of motion for <h_ is more complicated; at early times it can be expanded in powers of r/r[/ <C 1, 
where tjj is the age of the universe in conformal time. In what follows, we shall consider the leading 
order in such an expansion. 

Consider first the case of a radiation dominated Universe. 

For 4>_|_ we have the very same equation of GR, i.e. 


4 k'^ 

+ - 4>' + — $+ + O 
r 3 


m' 


■Ml 


pi 


= 0 . 


(37) 











For <I>_, outside the horizon where fcr <C 1 we get 


® $'_ + A _ 

r 


(k o? — /?^) 




( 38 ) 


and for the modes inside the horizon, namely /cr ^ 1, we have 


7 97 

T 11 


- 


(k — /5^) 
(k + /32) 




(39) 


Both scalars show no instability and their non-decaying modes are constant outside the horizon. In¬ 
side the horizon instead they both oscillate with a different decaying amplitude and frequency. The 
corresponding energy density perturbation is at the leading order 


P 


Iko^ 

3 (k -|- /J^) 


(3 + T^) $+ , 


(40) 


and therefore behaves like in GR. Outside the horizon it turns out to be constant, instead inside the 
horizon it has a fixed amplitude oscillating behaviour. In the case of matter dominated Universe both 
the fields <k+ and d>_ obey the same equation of GR at the leading order, namely 


4> 


+/- 




+ /- 


+ 0 ( — ) =0 
ru 


(41) 


It is worth to stress that the gauge-invariant density perturbation ^p^ijp is the real observable quantity 
and is given by 


P 


9 

Ka 

6 (ko^ -|- /I^) 


(12 +r^) $+. 


(42) 


Again the behaviour of bp^^jp is similar to GR: outside the horizon it is frozen and inside the horizon it 
grows like leading to the formation of structures. 

In order to study the other Bardeen potentials, it is useful to introduce also the combinations 


T+ = + 


K 


^ 2 , 


In general it turns out that 


T. = Ti - 


K 


4'2 • 


(43) 


T+ + 4>+ = 0 , 




rc p (k + /3^) ^ 


(44) 


It is interesting to note that in bigravity the double diff invariance is broken down to the diagonal diff 
invariance when an interaction between the two metric is introduced. The + sector is protected by 
diagonal diffs and indeed in the tensor sector corresponds to massless spin 2 modes, while the — gives 
rise to the massive ones, see also [39] . 

Remarkably, the instability that was present in the field 4>2 for modes inside the horizon with matter 
minimally coupled to found in [20], is not present. For such a behaviour the new matter coupling is 
instrumental, indeed the ^ ~ 1 regime for the background emerges and the very fast gradient instabilities 
disappear. However, as we will see soon, instabilities in the tensor and vector sectors emerge. 
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4.2 Tensor Sector 


The first hint of the problematic behaviour of fluctuations manifests itself when studying the tensor 
fluctuations. For this sector, the final expression for the perturbed Einstein equations is rather simple 


h-wpd(\ {hjfj - h^Jj) = 0 , 


tt' 


(45) 


+ 


2 |« + |)-^ 




2 ij 


+ {hlT,-h^)=0, 


(46) 


where /i, p and d are defined in Appendix B.l Also in this sector, it is useful to introduce the following 
combination of fields 

(3^ , ^ 


/i_(_ — hi -\- 


r ^2 


h- = hi — 


r ^2 , 


(47) 


with the indices and the symbol understood. Again it is h^ that is relevant for observed gravitational 
waves, since it is the combination appearing inside G^u- 

In particular, for the physical background solution, we get for general w 


h''+ 2 -Hh'+k^ /i+ = 0; 


(48) 


h'L + 2'Hh'_ + 


+ /3^) 


m 


fi w (^KO^ + p 


/32 


—a^ [ko^ — 0^) 


w? fi w[tia^ + l3P‘Yp 






h. = 0. 


h- 


(49) 


The combination /i_|_ behaves as in GR and represents spin 2 masseless modes protected by diagonal diff 
invariance. Instead the combination h- features a contribution to its mass proportional to w (see the 


parenthesis in the first line of eq (49)). This contribution is entirely due to the new matter coupling. If 


w is positive, the mode h- is tachyonic. If w is negative, this mode acquires a positive mass squared. 


Let us discuss some physical implications of the evolution equation (49) for the mode h-. During 


matter domination, u; = 0 and the new contributions to the effective mass to the mode /i_ vanish: there 
are no instabilities associated with the new coupling of gravity with matter. 

In the radiation dominated era, w = 1/3: the mode h- acquires a tachyonic mass, and an instability 
is expected. Indeed, radiation dominated super-horizon solutions for and h-, in terms of the scale 
factor a (X T, are 


h+ 

h_ 


a J 


K a 




{k o? + f30 


a 


VS+l 4 / 5-1 

Qj ^4“ G4 (x 2 


(50) 

(51) 


where Ci, C 2 , C 3 , are integration constants. The quantities C'i,C 2 control the healthy evolution of 


the mode /i+, and contribute also to the mode h- due to the source term in the second line of eq (49) 


The integration constants C 2 , C 3 correspond to decaying modes, that can be neglected, while controls 
a growing mode. As a consequence, this system features a mild power-like instability at superhorizon 
scales, but only during radiation dominated era. 

Hence, in order to ensure that this set-up is under perturbative control, we need to impose that the 
amplitude of /i_ does not exceed unity during radiation domination. Let us choose units in which the 
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scale factor after the reheating period, the beginning of radiation domination, is art = 1- Focus on the 
contribution of the growing mode only, and take into account the evolution of the scale factor during 
radiation domination. We find that at the end of the radiation period, at matter-radiation equality, we 
have to satisfy the inequality 


Ci {cLeq) 


Vs-i 

2 


C 4 



< 1 


(52) 


where the temperature of transition from radiation to matter domination is T^q — 3 eV, while the 
reheating temperature when the radiation era starts is model dependent, and depends on the reheating 
temperature after inflation. Taking a representative value Trh = 108 GeV, we find an upper bound for 
the integration constant 


Ci < 10 


-10 


(53) 


Interestingly, such small values for the quantity C 4 , can be motivated by the inflationary phase that 
precedes radiation domination. Indeed, it is inflation that sets the initial condition for the amplitude of 
tensor fluctuations, that evolve in the radiation dominated era. We sketch here an argument to explain 
this fact. Inflation is a phase of quasi-de Sitter expansion, where the parameter w ~ —1. During this 


phase, the mode h- acquires a positive mass squared in its evolution equation (49) due to the new 


coupling. The solution of the coupled system of equations (48), (49) is then (the scale factor scales as 
a ~ 


h+ 

h- 


D 2 


iyK o? — /3^) 


Ka- 


+ /32) 


Di- 


a3 


-|- ZI3 o 2 


sm 


'^/3 


In a ) -|- D 4 a 2 cos 


^%/3 


Ino 


(54) 

(55) 


with Di... D 4 integration constants. The quantities Di, D 2 control the evolution of the mode /i+; 
neglecting the decaying mode D 2 , the constant mode /i+ matches continuously between the inflationary 
and radiation dominated era, hence we set Di = Ci. The quantities D 3 , D 4 control the specific properties 
of the mode /i_, so they determine the initial conditions for the mode C 4 at the beginning of radiation 
era, after the end of inflation. Notice that ZI 3 , D 4 are decaying modes. We can write the matching 


relation between solutions (51) and (55) 


L >4 ~ C 4 , 


(56) 


where we make the crude assumption of a sudden transition from inflation to radiation domination, 
so that the end of inflation occurs at a^h = I- Assuming that during inflation the size of /i_ remains 


bounded, and that there are no cancelations among the different terms in eq (55), at the beginning of 
inflation we have the condition 

D 4 a-'/'<I. (57) 


Notice that in our units in which arh = I, the value of ain, the scale factor at the beginning of inflation. 


is very small. So in these units it is natural to choose a very small value for D 4 to satisfy condition (57). 
Indeed, saturating the previous equality, and using the relation arhldin = with W/ the e-fold 

number (that we take for definiteness N(.f = 60), recalling that in our units arh = I, we find 


D 4 = D4a2‘^ = ^ jQ-ag 


7)4 ~ (74 ~ 10 


-39 


(58) 


that comfortably satisfies inequality (53). 


II 













Of course the previous arguments are based on linear perturbation theory. Non-linear effects couple 
different modes: higher order fluctuations in the scalar sector are expected to feed the initial value 
of the tensor amplitude at the beginning of radiation era. On the other hand, in estimating their 
effects, one would have to take in careful consideration the coupling factors between scalars and /i_, 
that can be suppressed by factors of the graviton mass. So it is possible that even taking into account 
non-linearities, the growth of tensor fluctuations can still be maintained under control: it would be 
interesting to investigate more in detail this topic. 


4.3 Vector Sector 


The dynamics of fluctuations in the vector sector is unfortunately much more problematic. Using the 
conservation of the matter EMT with respect to the effective metric (equation ( |78[ ) in Appendix 
[a] ), the velocity perturbations 6vi can be completely solved. From equations (88 ](^ ) in Appendix |B.2[ one 
can show that all vectors can be expressed in terms of a single combination V 12 that satisfies a second 
order equation; thus, only a single transverse vector propagates. As for the tensors, the qualitative 
features of vector dynamics depends on the equation of state of the fluid constituting the matter EMT. 
When the fluid equation of state p = w p is such that w < 0, no instabilities are found and the system 
is healthy. Instead, we find serious gradient instabilities when tc > 0. In particular, in the case of a 
radiation dominated universe we have that 6v = 6vo, where 6vo is an arbitrary constant, and at the 
leading order in t/t (7 


Ui = - 


U2 = - 


Vi"o + 




Svo - 


5/32 


(k + j3‘^) -|- 5) 

5 K 

/c^ ^ ' {k a'^ + (3'^) {k‘^ + 5) 
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Svq + 


K 2 ; 

K 2 ; 


r (A:^ -|- 5) 


, (A:^ -I- 5) 

V(2- ^ ^ Vl2 = 0. 


5r^ 


(59) 

(60) 
(61) 


So we find the ‘wrong sign’ in front of the gradient term for the propagating mode V 12 , leading to a 
gradient instability in the sub-horizon limit kr ^ 1. In this regime, we have an exponential growth for 


V 12 , and the leading contribution to the solution of (61) is 


V12 oc e \/5 


(62) 


Also when selecting tuned initial conditions, such exponential instability is so severe that even a small 
initial amplitude of vector generated by non-linear effects (as the ones mentioned at the end of the 
previous section) will be rapidly amplified to a level that drives perturbation theory out of control. 
Notice that once more in the -|- sector no instability is present. Indeed, the combination 


U+ = Ui + 


/3' 


K 


V2 = - 


k? 


1 + 1 ^vo , 


K a‘- 


(63) 


represents a decreasing mode. Thus instabilities are present only in the — sector. To conclude, the 
dynamics of vector fluctuations ruin the cosmology of bigravity doubly coupled to matter. Possible 
ways-out, to be investigated in the future, could be to modify our ansatz for the EMT. Indeed in this 
work we mainly focussed on an EMT of perfect fluid form: it might be that other choices of matter 
content can lead to better behaved system. As an example, one could further include vector degrees 


‘‘Spatial indices for the vectors are understood. 


12 
















of freedom, also non-minimally coupled to gravity (as recently explored in massive gravity or related 
systems in @01 mills]), and study cosmological configurations in the scenario of doubly coupled matter. 
Acting as additional source, they might be able to fix the problems we found with vector fluctuations. 
We leave this interesting issues to a future investigation. 


5 Conclusions 


Massive gravity has been the subject of an extensive investigation. The main phenomenological mo¬ 
tivation is to explain the present acceleration of the universe. However devising a satisfactory model 
is not an easy task. The simplest ghost free Lorentz invariant version of massive gravity [35| with an 
auxiliary non-dynamical metric has no flat homogenous solution [7j and only Lorentz breaking mod¬ 
els [la El EH ESI ET] can support such configurations m- Sticking with Lorentz invariant models, to 
overcome the limitations of a non-dynamical auxiliary metric, one can promote it to a dynamical one 
in the contest of bigravity. As far as homogeneous configurations are concerned, in bigravity theories 
the situation drastically improves with respect to massive gravity, and one finds branches of flat FRW 


solutions. Leaving aside the ones with a curvature singularity at late time with c < 0 (see eq. (20)) - as 
and more recently in iiaEaEs] 


discussed in |T8| and more recently in - regular background solutions unfortunately features 

an exponential instability in the scalar sector at early time, quickly invalidating perturbation theory 
already during radiation domination [50]. Things do not change when both metrics are coupled with 
two different matter sectors [21]. In the presence of two metrics there is a certain degree of ambiguity 
on how to couple matter with gravity |24j and actually one can consider a sort of democratic coupling 
of the two metrics with matter, see eq. ([^. Though the new coupling reintroduce the Boulware-Deser 
ghost m. one can argue that its mass is above the cutoff and does not affect the low energy physics in 
the spirit of effective field theory. One of the interesting features of the new coupling is that it gives rise 


to an effective background dependent soft mass, see for instance (35) for the case of an FRW background. 
As a result even taking the “hard” mass m in the deforming potential to zero, we still have a massive 
gravity theory thanks to the environmental soft mass, that for an FRW background is proportional to the 
Hubble parameter. This is the feature that opens up a new dynamical regime compared with minimally 
coupled case such that ^ = u/a stays constant at early time and then flows toward a de Sitter attractor 
responsible for the present acceleration phase, see figure 1. The next step is to study the behaviour of 
cosmological perturbations. Things get better in the scalar sector where no instability is found, how¬ 
ever troubles develop in the tensor and specially in the vector sector. Among the tensor modes the -|- 


combination (see eq.(47)) which is protected by diagonal diffs has the same dynamics as in GR, while 
the second independent tensor mode develops a power-law growth until the matter domination. Such 
a growth is naturally counterbalanced by a sufficient low primordial production during inflation. In 
the vector sector the situation is more worrisome: while again the -|- combination of vector fields, see 


eq.(63), has only a decreasing mode, the — combination shows an exponential instability when w > 0 


at subhorizon scales. The main effect is to loose theoretical control of the theory at the perturbative 
level already during the radiation domination. On the other hand, we stress that in all cases studied the 
anomalous growth is only present in sectors that do not couple directly with observed matter; moreover, 
we speculated on possible extension of the bigravity theory under consideration, that might be able to 
cure such instability problems. 


®The term refers to the existence (non existence) of Lorentz global symmetry in the gravitational sector in the unitary 
gauge. 
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A Perturbed geometry 


fJ.U 1 


Let us now consider the perturbations of the FRW background (20) 

d/iu ~ dfiu T ® ) f/iu ~ ffiu T ^ 

parametrized as follows 

^100 = —2Ai , /i 2 00 = — 2 c^ A 2 , 

^l/20i = Ci/ 2 i — diBii 2 , d^Vi/ 2 i = d^Ci/ 2 i = = S^^h^'^i/ 2 ij = 0. 

hi/2ij = h^'^ii2ij + diVif2j + 9jVi/2i + 2didjEii2 + 2 5ij F 1/2 ■ 

Spatial indices are raised/lowered using the spatial flat metric. 

Under a gange transformation generated by the metric pertnrbation transforms 

Shi^u = a~'^ {Cdagfiu + gau d^C + g^la duC) , 

5h2^iv = [Cda Ifiu + fau d^C + Ua O^C) . 
and for the corresponding components 

SA, = nc^ + C^', 5B, = c^-c', 6E, = c, SF, = nc°-, 

6A2 = n0C^ + C^', <5i?2 = c2 c°-C', se2 = C, -5F2 = h.c°; 

<^Ci/2i = Ct ) <5V’i/2i = Cr ) Sh^'^ii2ij = 0; 


(64) 


(65) 


( 66 ) 


(67) 


where 


C = CT + ^^C, C = a-^9,C 

{coj)' d 

77/3 = 7—r = —I- 77a; • 

(cwj c 


( 68 ) 


In the scalar sector we have 8 fields and two independent gauge transformations, as a result we can form 
6 independent gauge invariant scalar combinations that we chose to be 


Ti — Ai — 77 Hi — H) 


T 2 = A 2 + c M - - 77, 


:: ‘='2 

^2 - 2 

& 


$1 = Fl - 77 “1 , $2 = 7^2 - 77a; ^ 


(69) 


,2 ’ 


£ = Ei-E2, Bi = B2- c^Bi + (1 - c^) E[ , 


where H 1/2 = 77 i /2 + By 2 - following additional gauge invariant fields will be useful to write in a 

compact form the perturbed Einstein equations 


-Hi = F 2 - Fi + (77 - 77a;) Hi, J-2 = F 2 - Fi + (77 - U^) Ha/c^ , 

B2 = B2-C^B^ + {1-^)E’2, 

Ai = c(A 2 — Ai) + [c (77 — 77a;) — c'] Hi, 

A 2 = c{A 2 — Ai) + [c (77 — 77a;) — c^] ^ 2 / 0 ^. 


(70) 
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The fields in (70) can be of course expressed in terms of the ones in (69). In the matter sector, since 
matter is minimally coupled to the effective metric we define the gauge invariant perturbed pressure 
and density in the following way 


where 


and 


^Pgi — Sp 


'-^eSP 


'aa + Pujc) 


2 ’ 


6p„i = 6p- 


-eSP 


(a a + /3 a; c) 


2 ’ 


Seff = BeS + E' 


eff 1 


BeS = Bi -\ -- {c Bi + B 2 ) + /3^ a;^ B 2 , 

1 + c 

E'eS = {oi a + j3 uj) (ya a E'i + j3 u] E'^ . 


(71) 

(72) 


(73) 

(74) 


For matter, together with pressure and density perturbation, there is also the perturbed 4-velocity 
that consists of a scalar part v and a vector part 5zi 


u^^ = u^^ + 5u^, u^u'^Gf,y = -l 


(75) 


The corresponding gauge invariant quantities are defined as 


Us = V + 


E' 


eff 


where 


{aa + j3 uj)^ ’ 

2 a(3 au 


Svi = Szi + 


^eSi 


CeSi = a^ Cli + 


1 + c 


(a a + /3 oj)^ ’ 
{cCli -\- C 2 i) + 0^ o;^ C 2 i ■ 


(76) 


(77) 


The conservation of the matter EMT with respect to the effective metric leads to the following 
differential relation for vector matter perturbations 


Sv'i — - - - -r 

{aa + /3u]) 


[{3 w — 1) {aaB + f3LdBuj) 


(3 uj {(3 uj o' + aa[d — {c — 1) {B — Buj)]} 

{aa + I3uj c) 

In the vector sector we have 4 fields and 1 gauge transformation; thus, we can form 3 independent gauge 
invariant vector perturbations 


dvi = 0 . (78) 


El 


= Ci /o,-V( 


1/2 i — ‘-'1/2 i 


l/2i 


Xi — I'l i ^2 i 


(79) 
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B Perturbed Einstein Equations 
B.l Scalars 


Using the definitions of the previous section we have for scalar perturbations of g 


2Ad>i + en {^in - ^[) + f 2 + pa (3 e] (3-Fi - ^8) 

Hip' 


aa^y^ 

L 




2 ^ 1 % - 2$; + ^ [rr? /s - y^ p] Bi + 


{1 + w) pay^ y 2 

Vc 


Us H- £ 


= 0 ; 


where 

and 


{didj - 5ijA) [a^ [firn^ - w pd^) £ - - ^i] 

+Sij [2a^ (m^ fi — w pd^) 3Fi + {m^ f 2 — w pa/3 y'^ Ai 


+2^1 + 2%') -2 - 2(2 

ao?y‘^ycW [ , 13ip' f „ 

= — ^{yiBi-y£} 

L Pc J 


/i = C PC (3 03 + 02 (c+ 1)) +oi] , /2 = c (hosC^ + 4o2C + ai) , 


y = {a + /3i) 
2 a 


yi = 

ys = 


+ I3i, 


1 + c 

al3y‘^i 
(1 + c) Vc 


yc = (a + /3c0, 
2l3ic 

2/2 = 0 + 


p = 




d = al3yyc, 


1 + c ’ 

[c y + re (a (1 + 2 c) + /? c (2 + c) ^)] . 


(80) 

(81) 


(82) 


(83) 


(84) 
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For the metric / we have 




2 c2 A^>2 + (^2^0; - ^ 2 ) + ^ (m2/2 + p « /? y' 0 {^£ ” 3 J 2 ) 


d' (3 y'^ 

L' 


^Psi + -^ {y2B2-c^y£') 


2 c (^ 2 ^^-^ 2 ) - 


(1 _|_ c) 


(m ^/2 - yz p) B 2 


+ 


d d {l + w) pPy^ yi 
K^Vc 


u,--£:'] = 0 : 

y 


-c {didj - 6ijA) 


+5, 




2d c 


(/i rd -wpdd) £ + c{^2 + 'I' 2 ) 
d c 


{rdfi -wpdi) T 2 H-^ {rdf 2 -wpal3y^^) A 2 

Ki, 


+2 ^2 ^ 2^' - 2 - U 2 - 2$" + 2 - - 2^^ U.' + 2^^ 'F' 


c 


2 M2 kC 


wf3y‘^yc6ij \ 6pgi + 5^(y2^2-c^y£:') 

^ Vc 


(85) 


( 86 ) 


(87) 


B.2 Vectors 

In the vector sector the perturbed Einstein equations are 

^dii 1 / 2f , -c \ (1 + w')y«y^y2 . „ 

A + f«!/4)x.-- Sv, = o 

{rd h -wpd^) %Vi 2 ,) = 0 ; 


AF 2 * ^ 1 


2dc ndi^dc) 


f 2 f , -c \ {l + w)pPy^yi 

[m h + P^ydXi -?- 0Vi = 0; 

I^^Vc 


%^ 2 j) + 


2{n+^]-- 


2 

d{iV2j) + ^{rd fi-wpd^) 9(iVi2j) = 0 
Kl; 


where 

and 


Vl2i = Vli-V2i, Ei2i = Fli-E2*, 


y4 = {2yc + ac + 13^ + wyc) ■ 

1 + c 


Notice that Fi 2 i = Xi ~ 1712 i- 


( 88 ) 

(89) 

(90) 

(91) 

(92) 

(93) 
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